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Abstract:  In  order  to  transmit  energy  in  one 
direction,  sonar  flextensional  transducers  are 
combined  into  arrays  of  elements  that  are  spaced 
a  1/4  wavelength  apart.  The  directionality  (front- 
to-back  pressure  ratio)  is  a  modest  of  6  dB.  A 
single  projector  that  is  1/3  wavelength  in  size 
and  capable  of  developing  unidirectional  beams 
with  front-to-back  ratios  greater  than  20  dB  that 
is  independent  of  sound  speed  is  described  here. 
The  directional  Class  VII  “dogbone” 
flextensional  is  a  modified  version  of  an 
elliptical  shaped  Class  IV  flextensional  having  a 
concave  shell  rather  than  a  convex  one.  The 
piezoelectric  ceramic  stack  is  a  tri-laminar  bar 
with  two  active  sections  separated  by  an  inactive 
center  section.  By  driving  both  sections  of  the 
stack  in-phase,  the  shell  is  driven  into 
omnidirectional  radiation  pattern.  Driving  each 
stack  section  1 80  degrees  out-phase  causes  a 
bending  mode  in  the  stack  resulting  in  a  dipole 
radiation  pattern.  By  driving  both  sections  using 
complex  coefficients  determined  from  the 
omnidirectional  and  dipole  mode  patterns,  a 
cardioid  directional  radiation  pattern  is 
developed.  Cardioid  patterns  are  developed  over 
an  octave  frequency  band  with  a  front-to-back 
pressure  ratio  of  more  than  50  dB.  COMSOL 
Finite  Element  Analysis  (FEA)  code  with 
Acoustics  Module  is  used  to  predict  in- water 
electroacoustic  performance  and  compared  with 
an  experimental  data. 
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1.  Introduction 

Sonar  transducers  are  used  for  transmitting  and 
receiving  acoustic  energy  for  the  purpose  of 
detection  and  location  of  underwater  objects.  In 
figure  1,  a  ship  is  towing  a  variable  depth  sonar 
that  contains  sound  projectors  used  to  transmit 
the  acoustic  energy  and  a  tethered  line  array  of 
hydrophones  used  to  receive  the  reflected  sound 
waves  from  an  underwater  object1.  These  sound 
projectors  are  usually  of  a  Class  IV  elliptical 
type  flextensional  sonar  transducer  that  has 
proven  to  be  rugged  high  power  compact 
omnidirectional  projectors2.  In  order  to  direct 
energy  in  a  direction  to  minimize  the 


reverberation  levels,  avoid  overloading  the 
receiver  line  array  preamplifiers,  and  removing 
the  left-right  ambiguity  problem  caused  by 
symmetrical  or  Omnidirectional  beams;  these 
sound  projectors  are  combined  into  two  line 
arrays  or  two  planar  arrays  of  projector  elements 
that  are  several  wavelengths  long  and  spaced  a 
1/4  wavelength  apart,  as  shown  in  figure  2. 
Theoretically,  two  point  sources  spaced  a  1/4 
wavelength  apart  from  one  another  with  one  of 
the  sources  driven  90  degrees  out  of  phase  with 
the  other  source  will  produce  a  cardioid  type 
beam  pattern  with  a  null  in  the  back  that  has  a  3- 
dB  down  beamwidth  of  180  degrees3,  as  shown 
in  figure  3a.  In  practice  this  is  not  the  case.  Two 
Class  IV  flextensional  transducers  spaced  1/4 
wavelength  apart  from  one  another  with  one  of 
the  sources  driven  90  degrees  out  of  phase  will 
produce  a  front-to-back  ratio  of  only  6  dB,  as 
shown  in  figure  3b. 

We  will  now  describe  the  Directional  Class  VII 
Dogbone  Flextensional  Sonar  Transducer 
projector  which  generates  directional  cardioid 
beams  from  a  single  transducer.  The  transducer 
is  1/3  wavelength  in  size  with  a  front-to-back 
ratio  greater  than  20  dB  and  independent  of 
sound  speed.  This  transducer  design  is  capable 
of  generating  an  omnidirectional  radiation  beam 
pattern  and  a  dipole  radiation  beam  pattern 
which  uses  superposition  to  create  a  cardioid 
unidirectional  beam.  The  synthesis  is  described 
in  figure  4,  where  the  omnidirectional  source  is 
represented  as  a  positive  pressure  radiator  and 
the  dipole  source  is  represented  by  positive  and 
negative  pressure  radiators,  which  when 
combined  form  a  cardioid  pattern.  Thus,  the 
acoustic  pressure  output  is  doubled  in  one 
direction  and  nulled  in  the  other  direction.  This 
type  of  pattern  generates  a  directivity  index  (DI) 
of  4.8  dB  and  has  a  3-dB  down  beamwidth  of 
130°.  We  are  proposing  that  a  single  line  array 
of  these  Directional  Class  VII  dogbone 
flextensional  transducers  could  replace  a  dual 
line  array  of  Class  IV  flextensional  transducers 
which  will  reduce  weight,  cost  and  increase 
front-to-back  ratio,  as  shown  in  figure  5.  We 
will  show  that  cardioid  type  patterns  can  be 
developed  over  an  octave  frequency  band  with  a 
predicted  front  (left)  to  back  (right)  pressure  ratio 
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of  more  than  50  dB  as  compared  with  6  dB  for 
the  1/4  wavelength  spaced  apart  projectors. 

2.  Description 

This  projector  merges  Directional  Flextensional 
Class  IV  Transducer  Technology4,5  with  that  of  a 
Class  VII  Dogbone  Flextensional  Transducer 
Technology6.  We  refer  to  this  as  the  Directional 
Class  VII  Dogbone  Flextensional  Sonar 
Transducer,  which  produces  enhanced  motion  on 
one  surface  and  canceled  motion  on  the  second 
surface.  The  Class  VII  flextensional  Sonar 
Transducer  is  referred  to  as  a  “dogbone” 
flextensional  because  of  its  similar  shape  to  that 
which  was  originally  described  by  Merchant  in 
U.S.  Patent  3,258,738  (28  June  1966).  This 
Class  VII  flextensional  is  a  modified  Class  IV 
shape  having  concave  staves  or  beams  radiating 
surfaces  rather  than  convex  staves  or  beams 
radiating  surfaces  (for  comparison  the  two  are 
shown  in  figure  6).  This  type  of  transducer  has 
the  advantages  that  all  the  surfaces  of  the  shell 
radiate  in  phase  and  a  positive  compressive  stress 
is  maintained  on  the  ceramic  stack  as  the 
operation  depth  increases.  The  shell  can  support 
thinner  walls,  since  hydrostatic  pressures  do  not 
have  to  be  over  come  in  order  to  keep  the 
piezoelectric  ceramic  stack  under  compression  as 
in  a  conventional  Class  IV  transducer. 

The  piezoelectric  ceramic  stack  is  a  trilaminar 
bender  bar  with  two  active  sections  separated  by 
an  inactive  section  which  is  inserted  along  the 
major  axis  of  the  shell.  Driving  both  sides  of  the 
stack  in-phase  allows  the  shell  to  be  driven  into 
its  conventional  mode  creating  an 
omnidirectional  radiation  pattern.  By  driving 
each  of  the  stack  sides  180  degrees  out-phase 
with  one  another  causes  a  bending  mode  in  the 
stack  creating  a  dipole  type  (cosine)  radiation 
pattern.  By  combining  these  two  modes  through 
superposition  generates  a  cardioid  type  radiation 
pattern.  The  uniqueness  of  this  transducer  is  its 
small  size  (less  than  a  third  of  a  wavelength) 
producing  directional  beam  patterns  at  high 
acoustic  output  power  from  a  single  element, 
unlike  1/4  wavelength  dual  line  arrays. 
COMSOL-Multiphysics  finite  element  code 
Acoustics  Module7  with  a  Perfectly  Matched 
Layer  ( PML )  as  a  boundary  condition  is  used  to 
predict  in- water  electroacoustic  performance 
(Impedance,  TVR  and  Beam  Patterns)  of  a  Class 
VII  dogbone  transducer  and  compared  with  an 
omnidirectional  3  kHz  Class  VII  flextensional 
transducer  experimental  prototype  unit. 


COMSOL  is  then  used  to  predict  the  far  field 
complex  pressures  of  the  omnidirectional  and 
dipole  modes  which  are  used  to  determine  the 
complex  drive  coefficients  used  to  drive  the 
transducer  into  the  directional  mode. 

3.  Model 

Measured  results  from  an  omnidirectional  3  kHz 
Class  VII  flextensional  transducer  experimental 
prototype  unit  will  be  used  to  validate  a  2-D 
FEA  model  before  developing  a  directional  Class 
VII  flextensional  transducer  model.  The 
experimental  transducer  has  two  shells  and 
stacks  electrically  wired  in  parallel  and  two  end 
plates  separated  by  four  long  bolts.  It  is  then 
encapsulated  in  a  water  tight  polyurethane  boot. 
Photos  of  the  transducer  Shell  and  Ceramic 
Stack  are  shown  in  figure  7a  and  encapsulated  in 
figure  7b.  Figure  8  shows  the  dimensions  of  the 
transducer  model  which  is  composed  of  a 
piezoelectric  drive  stack  and  vibrating  shell.  The 
stainless  steel  shell  has  thick  rigid  ends  which 
the  piezoelectric  drive  stack  pushes  against 
which  in  turn  flexes  the  curved  shell  staves  (or 
Beams).  The  staves  are  4.32  mm  thick  and  have 
an  inner  and  outer  radius  of  85.85  mm  and  90.17 
mm  respectively.  The  piezoelectric  ceramic 
drive  stack  is  composed  of  30  PZT-8  U.S.  Navy 
type  III  ceramic  3  3 -mode  plates  that  are  each 
2.54mm  thick,  17.78  mm  wide  and  40.82  mm 
deep,  two  of  which  are  unpoled  and  used  for 
insulators  on  the  ends.  Figure  9  shows  the  2-D 
FEA  Stack  Subdomain  Settings  and  the 
Boundary  Condition  Settings  for  Piezoelectric 
Plane-Strain  application  FEA  model,  Plane- 
Stress  FEA  modeling  for  Piezoelectric 
application  is  not  available.  The  glue  joints  and 
brass  foil  electrodes  that  are  between  each  plate 
that  are  used  to  consolidate  and  wire  the  stack 
are  not  used  here  in  the  model.  Instead,  the 

e 

piezoelectric  material  compliance  matrix  s33  was 

modified  from  13.5xl0'12  m2/N  to  21xl0'12  m2/N 
to  account  for  the  added  compliance  of  the  glue 
joints  and  foil  electrodes,  which  lowers  the  stack 
in-air  resonance  from  22  kHz  to  that  of  the 
measured  stack  in-air  resonance  of 
approximatlyl7.3  kHz  before  it  was  inserted  into 
the  shell. 

4.  2-D  Ceramic  Stack  Model 

As  shown  in  figure  9  Stack  Subdomain  Settings, 
the  ceramic  plates  are  labeled  PZT8U  and 
PZT8D.  This  is  to  account  for  the  electrical 


polarization  field  polarity  pointing  to  the 
direction  of  the  positive  electrode  surface  of  the 
ceramic  plate,  thus  the  coordinate  system  must 
be  rotated  1 80  degrees  to  one  another.  Figure  9c 
shows  the  electric  potential  and  field  arrows 
surface  plot  results  with  a  positive  rms  voltage 
potential  of  1  applied  to  the  boundary  settings 
listed  in  figure  9b.  COMSOL  does  not  allow 
direct  extraction  of  electrical  current  under  the 
piezoelectric  application  in  order  to  calculate  the 
electrical  input  impedance  (Z=V/I)  or  electrical 
input  admittance  (Y=I/V)  that  are  used  to 
characterize  sonar  transducers.  Thus,  by 
integrating  the  Normal  Current  Density 
(nJ_smppn)  or  Displacement  Current  Density 
(Jdy_smppn)  with  respect  to  the  ceramic  plates 
surface  area  the  inward  current  can  be  calculated. 
The  Normal  Current  Density  would  be  used 
when  there  are  glue  joints  between  plates  or 
having  an  uncommon  potential.  Displacement 
Current  Density  would  be  used  when  there  are 
no  joints  between  plates  and  have  a  common 
potential,  which  is  our  case  here.  The 
expressions  in  equations  (1)  and  (2)  are  the 
Displacement  Current  Density  assigned  to  each 
of  the  positive  potential  (+V)  surfaces  listed  in 
figure  9b.  The  Boundary  Condition  Settings  the 
“up”  and  “down”  in  the  equations  represents  the 
direction  of  the  current  flow  that  is  associated 
with  the  chosen  ceramic  plates  polarity  field  of 
the  PZT8U  and  PZT8D  plates. 


IU  =-(up(Jdy_smppn))*z 

(i) 

ID  =-(down(Jdy_smppn))*z 

(2) 

The  “z”  represents  the  depth  of  the  ceramic 
plates  40.82  mm  and  IU  is  the  assigned  the  name 
for  Equation  (1)  and  ID  is  the  assigned  the  name 
of  equation  (2),  once  these  expression  are 
integrated  under  COMSOL’ s  Integration 
Coupling  Variable>Boundary  Integration 
Variables  menu  Option.  Therefore,  by  summing 
the  magnitudes  of  both  IU  and  ID  is  the  total 
electrical  current  into  the  piezoelectric  stack. 

The  admittance  magnitude  is  given  by  the  sum  of 
the  magnitudes  of  both  IU  and  ID  divided  by  the 
applied  electrical  potential  (V=l),  therefore  the 
admittance  magnitude  is  given  by  equation  (3). 

|  Y|=I/V =abs(ID)+abs(IU)  (3) 

Figure  10  shows  the  piezoelectric  stack  in-air 
modeled  admittance  magnitude  response  of  the 
model  shown  in  figure  9  using  equation  (3) 
compared  with  a  measured  stack  that  was  used  in 


the  prototype  transducer  shown  in  figure  7. 
Admittance  is  a  complex  quantity  Y=G+jB, 
where  G  the  conductance  is  the  real  part  and  B 
the  susceptance  is  the  imaginary  part  of  the 
admittance,  they  are  given  in  equations  (4)  and 
(5)  for  V=l, 


G=  real(ID)-real(IU) 

(4) 

B=  imag(ID)-imag(IU) 

(5) 

The  IU  is  subtracted  from  ID  to  account  for  the 
opposite  direction  of  the  current  flow  that  is 
associated  with  the  chosen  ceramic  plate 
polarity,  the  order  of  subtraction  maybe  different 
for  other  models  due  to  boundary  condition 
choices.  The  capacitance  is  given  by  B/oo,  where 
co=27if  and  therefore  is  given  by  equation  (6) 

C=  (imag(ID)-imag(IU))/(2*pi*freq)  (6) 

where  “freq”  is  the  frequency  in  Hz.  The  low 
frequency  1  kHz  capacitance  of  the  stack 
predicted  from  the  model  was  64  nFd  verses  the 
measured  77  nFd.  The  models  relative  free 

T 

dielectric  constant  £r  is  1000  verses  that  of  the 
manufactures  data  sheet  EDO  EC-69  value  of 
1050. 

5.  2-D  Shell  and  Stack  Model 

The  structural  Subdomain  Settings  and  the 
Boundary  Condition  Settings  for  the  Shell  and 
Piezoelectric  Stack  are  shown  in  figure  1 1 .  The 
Shell  material  properties  are  of  stainless  steel  and 
the  piezoelectric  stack  is  the  same  shown  in 
figure  9.  The  shells  are  each  50.8  mm  in  height 
and  the  stacks  are  each  40.82  mm  in  height.  The 
plane-strain  application  in  COMSOL  Acoustics 
Module  allows  a  thickness  of  a  material  to  be 
entered  into  the  z-  axes,  which  otherwise  it 
would  be  a  unit  or  1  meter  thickness.  For  the 
Shell/Stack  model  0.1016  m  was  entered  for 
shell  thickness  and  0.08164  m  for  stack 
thickness,  1  meter  and  0.80354  meters  for  shell 
and  stack  respectively  would  also  work.  For  the 
proper  current  results  the  “z”  value  in  equations 
1  and  2  is  0.08164  m.  The  fluid  Subdomain 
Settings  and  the  Boundary  Condition  Settings  for 
pressure  acoustics  (acpr)  application  mode  are 
shown  in  figure  12.  The  fluid  loading  is  in  two 
parts  the  inside  water  boundary  and  the  outside 
boundary  the  Perfectly  Matched  Layer  ( PML ). 
The  water  boundary  contains  the  acoustic 
properties,  the  density  and  speed  of  sound  of  the 


(11) 


fluid  and  is  two  wavelengths  in  radius  at  3  kHz 
or  1  meter  and  is  considered  a  far- field  radiation 
condition,  this  is  where  the  acoustic  information 
is  contained.  The  PML  is  the  radiation  condition 
for  perfectly  absorbing  plane,  cylindrical  and 
spherical  waves  and  absorbs  incident  radiation 
without  producing  reflections,  here  the  model  is 
0.5  meters  or  one  wavelength  wide.  The  PML  is 
a  cylindrical  type  loading  for  this  2-D  geometry 
model,  thus  the  acoustic  pressure  will  decrease 
cylindrically  in  the  radial  and  z-axis  direction. 

At  the  interface  between  the  water  and  shell  an 
outward  pressure  -p  in  terms  of  force  per  unit 
area  is  applied  as  a  Load  around  the  surface  of 
the  shell  and  is  given  by  equations  (7)  and  (8).  A 
normal  acceleration  an  is  applied  around  the 
surface  of  the  shell  in  equation  (9);  this  boundary 
condition  is  used  to  couple  the  acoustic  domain 
to  structural  domain.  A  far  field  p^far  condition 
is  applied  around  the  outer  boundary  of  the  PML. 

Fx=  -p  *nx_smppn  (7) 

Fy--p  *ny_smppn  ( 8 ) 

an  -  nx_smppn  *u_tt_smppn + ny_smppn  *v_tt_smppn  (9  ) 

The  2-D  FEA  Fluid  and  Shell/Stack  Mesh  is 
shown  in  figure  13,  the  mesh  contains  16312 
triangular  elements.  The  typical  rule  of  thumb 
for  element  size  is  greater  then  five  elements  per 
wavelength7.  This  is  to  allow  higher  frequency 
calculations  and  as  can  be  seen  we  are  far  greater 
then  that.  COMSOL  application  scalar  variables 
pressure  reference  pref  should  be  change  to  IE-6 
Pa  this  is  the  acoustic  sound  pressure  reference 
level  for  water.  Also  both  the  excitation 
frequency  variables  freq_smppn  and  freq_acpr 
should  contain  the  expression  freq  which  will 
couple  the  acoustic  domain  and  structural 
domain  together  for  performing  a  frequency 
response  and  a  time  harmonic  analysis. 

COMSOL  acoustic  sound  pressure  level  (SPL), 
Lp  at  some  point  in  the  water  acoustic  pressure 
field  is  given  by  equation  (10) 

SPL  =10logw(pp*  /2Pref2)  (10) 

where  pref  is  the  acoustic  sound  pressure 
reference  level  for  water  1  uPa.  Here  the 
pressure  p  is  a  peak  value  since  the  2  in  the 
denominator  normalizes  pp  to  a  rms  value.  The 
more  common  form  in  terms  of  rms  pressure  is 


SPL  =  20  log  ,0  (/>„,„/  pnf) 

Transmit  voltage  response  (TVP)  is  an  important 
transducer  calibration  parameter  that  describes 
the  SPL  of  a  transducer  at  1  meter  distance  when 
driven  with  1  volt  rms  over  a  frequency  band  of 
interest  and  has  the  units  of  dB//l  uPafVrms  at 
lm.  Thus  to  determine  the  proper  TVR  level 
when  using  COMSOL’ s  SPL  function  a  voltage 
potential  of  V2  Volts  is  used  instead  of  1  Volt , 
since  a  peak  voltage  will  generate  a  peak 
pressure.  If  a  voltage  potential  of  1  Volt  was 
used,  3  dB  would  be  needed  to  be  added  to  the 
SPL.  As  discussed  earlier  the  2-D  model  has 
cylindrical  loading,  thus  both  the  acoustic 
pressure  will  decrease  cylindrically  in  the  radial 
and  z-axis  direction.  The  TVR  given  in  equation 
(12)  accounts  for  cylindrical  loading  effects  of 
both  the  radial  and  z-axis  directions 

TVR  =  SPL-10log10(l/h)-10log10(l/r)  (12) 

where  h  is  the  transducer  height  in  the  z- 
direction  and  r  is  the  radial  distance  pressure 
point.  The  transducer  is  0.1016  m  in  height  and 
radial  distance  pressure  is  1  meter.  The  TVR 
equation  is  now  TVR=SPL-9.93  dB. 

6.  In- Water  Results 

The  in- water  measured  and  modeled  transmit 
voltage  response  using  equation  (12)  is  shown  in 
figure  14.  The  resonance  at  3  kHz  is  the  first 
fundamental  flexural  mode  of  the  transducer 
with  a  design  of  operation  of  2.5  kHz  to  5  kHz. 
The  second  resonance  at  approximately  10  kHz 
is  the  breathing  (Hoop)  mode  of  the  Class  VII 
shell.  This  mode  happens  when  the 
circumference  of  the  ring  (beams)  are  equal  to 
the  wavelength  in  the  shell  material.  The  in¬ 
water  measured  and  modeled  conductance  G 
response  using  equation  (4)  and  susceptance  B 
response  using  equation  (5)  are  shown  in  figures 
15  and  16  respectively  for  voltage  potential  of 
1  Volt  applied  to  the  piezoelectric  stack.  The 
modeled  response  resonant  peaks  in  the 
conductance  response  are  lower  then  in  the 
measured.  This  is  a  result  the  cylindrical  loading 
which  will  overly  dampen  the  model  versus  that 
of  spherical  loading.  The  measured  and  modeled 
acoustic  radiation  beam  patterns  at  3  kHz,  5  kHz 
and  9.6  kHz  are  shown  in  figure  17.  The  acoustic 
beam  patterns  in  COMSOL  are  recorded  by 
plotting  the  SPL  verses  bearing  angle 


(atan2(y,x))  around  the  shell  at  the  radial 
distance  pressure  of  1  meter. 

7.  Directional  Dogbone  Model 

The  directional  dogbone  model  shown  in  figure 
18  is  a  slight  modification  of  the  model  shown  in 
figure  8.  Here  the  piezoelectric  ceramic  stack  is 
divided  into  two  active  sections,  a  Side- A  and  a 
Side-B  separated  by  an  inactive  section.  This 
will  allow  the  transducer  to  be  driven  into  the 
dipole  and  directional  modes.  This  is 
demonstrated  in  figure  19  by  applying  the  same 
voltage  potential  to  both  sides  of  the  stack.  This 
allows  the  shell  to  be  driven  into  its  conventional 
mode  where  both  beams  flex  in  and  out  with  the 
same  phase.  This  will  create  an  omnidirectional 
radiation  pattern,  see  figure  19a.  Applying 
opposite  voltage  potentials  to  each  of  the  stack 
sides  will  causes  a  bending  mode  in  the  stack 
which  will  cause  the  beams  to  flex  out  of  phase 
with  each  other  creating  a  dipole  type  (cosine) 
radiation  pattern,  see  figure  19b.  By  combining 
these  two  modes  through  superposition  by 
driving  the  stacks  with  complex  voltage  potential 
coefficients  will  minimize  motion  on  one  side 
and  enhanced  motion  on  the  other  side.  This  will 
generate  a  cardioid  type  radiation  pattern. 

8.  Complex  Voltage  Potential  Coefficients 

The  complex  voltage  potential  drive  coefficients 
for  the  transducer  to  generate  a  cardioid  type 
radiation  pattern  are  determined  from  acoustic 
SPL  amplitude  and  phase  of  each  of  the  modes  in 
the  far- field  and  then  combining  them.  The 
omnidirectional  mode  is  excited  by  applying  1 
volt  to  sides- A  ( Ea=l )  and  1  volt  to  side-B 
( Eb=l );  this  generates  a  far-field  complex  sound 
pressure  P0.  The  dipole  mode  is  excited  by 
applying  1  volt  to  sides-A  ( Ea=l )  and  -1  volt  to 
side-B  (. Eb=-1)\  this  generates  a  far-field 
complex  sound  pressure  Pd.  The  two  pressures 
should  always  be  at  the  same  angular  position. 
The  two-mode  synthesis  for  determining  the 
complex  voltage  potential  drive  coefficients  Ea 
and  Eb  in  algebraic  form  is  as  follows: 

Ea=(Pd+Pa)/Pd  =1  +  *  >  (13) 

Eb=(Pd-P0)/Pd=l~R  (14) 

where  R  =  Po/  p  .  For  enhanced  pressure  on 

side-A  the  omnidirectional  pressure  is  added  to 
the  dipole  pressure  as  shown  in  equation  (13) 


and  for  a  null  on  side-B  of  shell  the 
omnidirectional  pressure  is  subtracted  from  the 
dipole  pressure  as  shown  in  equation  (14).  Since 
each  of  the  modes  produces  different  SPL’s, 
equations  (13)  and  (14)  are  normalized  by  the 
dipole  mode  pressure  Pd.  Also,  since  these 
equations  contain  complex  pressures  with 
different  amplitude  and  phases,  the 
omnidirectional  pressure  phase  is  referenced  to 
the  dipole  pressure  phase.  The  ratio  of  equations 
(13)  and  (14)  is  given  by  equation  (15) 

Ea  /  Eb  =1  +  R/1-R  (15) 

which  normalizes  the  voltage  potential  of 
ceramic  stack  side-A  to  that  of  ceramic  stack 
side-B  in  terms  of  a  complex  quantity  or 
amplitude  and  phase  quantity.  Figure  20  shows 
the  a)  amplitude  pressure  and  b)  phase  plots 
angle  in  the  omnidirectional  and  dipole  modes 
that  were  used  to  calculated  the  complex  voltage 
potential  coefficients  at  3  kHz.  The  peak 
omnidirectional  pressure  is  Ppk  =  32. 73  Pa  and 
phase  is  <f>0=84.86  °,  therefore  the  rms  pressure  is 
Prms  =23.15  Pa  and  the  SPL  is  147.29  dB  using 
equation  (11)  and  the  TVR  is  137.36 
dB//luPa/V@lm  using  equation  (12).  The  peak 
dipole  pressure  is  Ppk  =  6.97  Pa  and  phase  is 
(j)d=- 144.36  °  therefore  the  rms  pressure  is  Prms  = 
4.93  Pa  and  the  SPL  is  133.86  dB  using  equation 
(11)  and  the  TVR  is  123.93  dB//luPa/V@ lm 
using  equation  (12).  These  TVR’s  were  then 
converted  back  to  pressures,  then  R  is  computed 
and  the  ratio  of  Ea/Eb  is  solved.  The  reasoning 
for  computing  the  TVR  is  this  is  the  most 
commonly  measured  quantity  on  underwater 
sonar  transducers.  The  complex  voltage  potential 
coefficients  for  2,  3  and  4  kHz  are  shown  in 
Table  I.  The  complex  values  were  applied 
instead  of  amplitude  and  phase,  side-A  voltage 
potential  is  Ea=(-0.721-j*0.244)*1.414  and  side- 
B  voltage  potential  is  Eb=l*1.414  at  3  kHz. 

Table  I.  Complex  Voltage  Potential  Coefficients 


for  2,  3  and  4  kHz. 
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The  modeled  pressure  surface  and  SPL  surface 
plots  in  the  omnidirectional,  dipole  and 
directional  modes  at  3  kHz  are  shown  in  figure 
21a  and  21b  respectively.  The  molded  beam 
pattern  omnidirectional,  dipole  and  directional 
modes  at  3  kHz  are  shown  in  figure  22a,  22b  and 
22c  respectively.  Using  these  coefficients  in 


Table  I  will  point  the  beam  to  side-B  direction 
instead  of  side-A  direction.  This  may  be  caused 
by  the  dipole  beam  patterns  lobe  on  side-B  180 
degrees  out-of  phase  with  side  -A  lobe  as  shown 
in  figure  4,  for  presentation  purposes  the  cardioid 
beams  were  reversed  to  point  to  side-A  direction 
in  figures  21  and  22.  The  cardioid  type  beam 
pattern  generated  in  figure  22c  has  a  67.6  dB 
front- to-back  ratio  (or  deep  back  null)  and  a  -3 
dB  beamwidth  of  142.5  degrees.  Cardioid  beam 
pattern  were  generated  from  1  kHz  to  6  kHz  in 
steps  of  100  Hz.  The  lowest  front- to-back  ratio 
is  56.5  dB  and  -3  dB  beamwidth  ranged  from 
131  degrees  at  1  kHz  to  164  degrees  at  4.5  kHz. 
The  TVR  omnidirectional,  dipole  and  directional 
modes  from  1  kHz  to  6  kHz  are  plotted  in  figure 
23  for  comparison.  The  back  null  side  of  the 
directional  mode  is  also  plotted  to  show  the 
relative  front- to-back  ratio. 

The  oceans  sound  speed  and  density  profiles8 
both  vary  with  water  depth.  Thus  it  would  be 
interesting  to  see  what  would  happen  to  the 
front-to-back  ratio  when  using  the  same 
coefficients  that  were  determined  with  the  sound 
speed  1500  m/s  and  density  1000  kg/m3  at  3 
kHz.  Figure  24  shows  the  front-to-back  ratio 
comparison  between  two  element  Class  IV 
flextensional  transducers  spaced  1/4  wavelength 
apart  shown  in  figure  3b  and  the  Class  VII 
Directional  Dogbone  flextensional  using  the  3 
kHz  coefficients  in  Table  I  for  different  sound 
speeds  while  keeping  the  density  constant.  The 
other  case  is  varying  the  sound  speed  while 
keeping  the  characteristic  impedance  (density 
times  sound  speed)  the  same  Rho-C  =  1.5  x  106 
Pa-m/s.  The  two  element  Class  IV  flextensionals 
had  a  constant  front-to-back  ratio  of  6  dB  from 
1300  to  1700  m/s,  while  the  dogbone  had  a 
greater  than  20  dB  front-to-back  ratio  for  both 
cases. 

9.  Conclusions 

COMSOL-Multiphysics  finite  element  code  and 
Acoustics  Module  with  a  PML  as  a  boundary 
condition  is  used  to  predict  in- water 
electroacoustic  performance  (Impedance,  TVR 
and  Beam  Patterns)  of  a  directional  Class  VII 
dogbone  transducer.  The  model  was  compared 
with  an  omnidirectional  3  kHz  Class  VII 
flextensional  transducer  experimental  prototype 
unit  with  good  agreement.  COMSOL  is  then 
used  to  predict  the  far  field  complex  pressures  of 
the  omnidirectional  and  dipole  modes  which  are 
then  used  to  determine  the  complex  drive 


coefficients  used  to  drive  the  transducer  into  the 
directional  mode  that  generated  cardioid  type 
beam  patterns.  We  showed  that  cardioid  type 
patterns  can  be  developed  over  an  octave 
frequency  band  with  a  predicted  front-to-back 
pressure  ratio  of  more  than  50  dB  as  compared 
with  6  dB  for  the  1/4  wavelength  spaced  apart 
projectors.  Thus  a  single  line  array  of  these 
Directional  Class  VII  dogbone  flextensional 
transducers  could  replace  a  dual  line  array  of 
Class  IV  flextensional  transducers,  which  will 
reduce  weight,  cost  and  increase  front-to-back 
ratio. 
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a)  Omni  b)  Dipole  c)  Cardioid 


Figure  1.  Sonar  uses  Transmitted  and  Reflected  Sound 
Waves  to  Locate  Underwater  Objects 


Figure  4.  Synthesis  of  a  Directional  Flextensional 
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F/B  <  15  dB 


b)  Directional 
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Figure  2.  Two  Line  Arrays  and  Planar  Arrays  of 
Projector  Elements  that  are  Several  Wavelengths  Long 
and  Spaced  a  1/4  Wavelength  Apart. 


Figure  3.  Modeled  Beam  Patterns  for  Two  sources 
spaced  a  1/4  wavelength  apart  with  one  of  the  sources 
driven  90  degrees  out  of  phase  with  the  other  sources, 
a)  for  Point  Source  Element  and  b)  Real  Case 
Flextensional  Elements. 


Figure  5.  Two  Line  Arrays  of  Projector  Elements  that 
are  Several  Wavelengths  Long  and  Spaced  a  1/4 
Wavelength  Apart  Replaced  by  One  Line  Array  of 
Directional  Dogbone  Transducers  Producing  a  Front- 
to-back  ratio  greater  then  20  dB. 


a)  Class  IV  Flextensional 


b)  Class  VII  Flextensional 


Figure  6.  Class  IV  Flextensional  Transducer  and  Class 
VII  “Dogbone”  Flextensional  Transducer 


a)  Shell  and  Ceramic  Stack 


b)  Encapsulated 


Figure  7.  Class  VII  “Dogbone”  Flextensional  Sonar 
Transducer,  Shell  and  Ceramic  Stack  and  b) 
Encapsulated 
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Figure  10.  Piezoelectric  Stack  In-air  Admittance 

Magnitude  Response  Measured  ( _ )  and  COMSOL 

FEA  Modeled  (oooo) 


Figure  8.  Class  VII  Dogbone  Flextensional 
Sonar  Transducer  Dimensions  in  mm. 


a)  Subdomain  Settings 


b)  Boundary  Settings 


a)  Subdomain  Settings 


Figure  9.  COMSOL  Stack  a)  Subdomain 
Settings  and  b)  Boundary  Settings  for  Piezo 
Plane  Strain  (smppn)  application  and  c)  Electric 
Potential  and  Field  Arrows  results 


Figure  11.  COMSOL  Shell  and  Stack  a) 
Subdomain  Settings  and  b)  Boundary  Settings 
For  Piezo  Plane  Strain  (smppn)  application 


Figure  12.  COMSOL  Fluid,  Shell  and  Stack  a) 
Subdomain  Settings  and  b)  Boundary  Settings  for 
Pressure  Acoustics  (acpr)  application 


16,000 


Figure  13.  COMSOL  2-D  FEA  Fluid  and  Shell/Stack 
Mesh 


Figure  14.  Transmit  Voltage  Response  Measured 
( _ )  and  COMSOL  Modeled  (oooo) 


Frequency  (Hz) 

Figure  16.  In- water  Susceptance  Response  Measured 
( _ )  and  COMSOL  Modeled  (oooo) 


Figure  17.  Beam  Patterns  in  Omnidirectional  Mode  a) 
Measured  and  b)  COMSOL  Modeled  at  3  kHz,  5  kHz 
and  9.6  kHz 


Figure  15.  In- water  Conductance  Response  Measured 
( _ )  and  COMSOL  FEA  Modeled  (oooo) 


Figure  18.  Directional  Dogbone  Flextensional 
Sonar  Transducer 
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Figure  19.  COMSOL  FEA  Piezoelectric  Ceramic 
Stack  Electric  Fields  and  Displacement  Shell  and 
Stack  Modes  for  Omnidirectional,  Dipole  and 
Directional  Modes  at  3  kHz 


Omni  Mode  Dipole  Mode 


A  A 

1 

z. 

1 

1' 

j  \ 

V  v 

a)  Peak  Total  Acoustic  Pressure 

J — -tZi—J — pJL- L 

j: 

2 

1  nt~ 

b)  Acoustic  Pressure  Phase  Angle 


Figure  20.  COMSOL  FEA  Molded  a)  Peak  Acoustic 
Pressure  and  b)  Acoustic  Pressure  Phase  Angle  in  the 
Omnidirectional  and  Dipole  Modes  at  3  kHz 


Figure  22.  COMSOL  FEA  Molded  Beam  Pattern  a) 
Omnidirectional,  b)  Dipole  and  c)  Directional  Modes 
at  3  kHz 
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Figure  23.  COMSOL  FEA  Modeled  Transmit  Voltage 
Response  Omnidirectional  Mode  (0000),  Dipole  Mode 

(AAAA),  Front  Side  Directional  Mode  ( _ )  and 

Back  Side  Directional  Mode  ( - ) 


Omni  Mode  Dipole  Mode  Directional  Mode 


a)  Pressure  Surface  Plots 


b)  SPL  Surface  Plots 


Figure  21.  COMSOL  FEA  Molded  a)  Pressure  and  b) 
SPL  Surface  Plots  in  the  Omnidirectional,  Dipole  and 
Directional  Modes  at  3  kHz 


Figure  24.  Front-to-back  Ratio  Comparison  Between 

Two  Lines  ( - )  and  Directional  Dogbone  ( _ ) 

and  Constant  Characteristic  Impedance  (0000)  using 
the  3  kHz  Coefficients  for  different  Sound  Speeds. 
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•  Sonar  transducers  are  electro-acoustic  devices  used  for  transmitting 
and  receiving  acoustic  energy  for  the  purpose  of  detection  and  location 
of  underwater  objects. 

•  In  order  to  transmit  energy  in  one  direction,  sonar  Class  IV  flextensional 
transducers  are  combined  into  arrays  of  elements  that  are  spaced  a  1/4 
wavelength  apart.  The  directionality  (front-to-back  pressure  ratio)  in 
practice  is  a  modest  6  dB  due  to  diffraction. 

•  A  new  class  of  transducer  the  Directional  “dogbone”  flextensional 
transducers  which  generates  cardioid  beams  could  replace  these  dual 
line  array.  This  will  reduce  weight,  cost  and  increase  front-to-back  ratio 
greater  than  20  dB. 

•  COMSOL  with  Acoustics  Module  is  used  to  predict  in-water  electro¬ 
acoustic  performance  and  is  compared  with  experimental  data. 
COMSOL  is  then  used  to  calculate  the  complex  drive  coefficients  used 
to  drive  the  transducer  into  the  directional  mode. 
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Two  Line  Arrays  and  Planar  Arrays  of  Projector  Elements 
that  are  Several  Wavelengths  Long  and  Spaced  a  1/4 

Wavelength  Apart. 
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Modeled  Beam  Patterns  for  Two  sources  spaced  a  1/4 
wavelength  apart  with  one  of  the  sources  driven  90  degrees 
out  of  phase  with  the  other  sources,  a)  for  Point  Source 
Element  and  b)  Real  Case  Flextensional  Elements. 
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An  other  way  to  generate  a  directional  cardioid  Beam 
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Two  Line  Arrays  of  Projector  Elements  that  are  Several 
Wavelengths  Long  and  Spaced  a  1/4  Wavelength  Apart 
Replaced  by  One  Line  Array  of  Directional  Dogbone 

Transducers 


a)  Conventional 
two  lines 
F/B  <  15  dB 


b)  Directional 
one  line 
(Replacement) 
F/B  >  20  dB 
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Class  IV  Flextensional  Transducer  and  -1 
Class  VII  “Dogbone”  Flextensional  Transducer 


a)  Class  IV  Flextensional 


b)  Class  VII  Flextensional 


NEWPORT 


Model  Verification 
Procedure 


Revolutionary  .  .  .  Releront  Results 


•  Validate  COMSOL  Model  with  Measured 
Results  of  a  standard  Omni  Mode  Dogbone 
Flextensional  Sonar  Transducer 

•  Then  predict  the  directional  operation  mode 
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NEWPORT 

Class  VII  “Dogbone”  Flextensional  Sonar  Transducer, 
a)  Shell  and  Ceramic  Stack  and  b)  Encapsulated 


a)  Shell  and  Ceramic  Stack 


b)  Encapsulated 


Page  10 


WARFARE  CENTERS 


NEWPORT 


Revo! ut  Umar  it  Research  .  . 


COMSOL  Stack  a)  Subdomain  Settings  and  b)  Boundary 
Settings  for  Piezo  Plane  Strain  (smppn)  application  and  c) 
Electric  Potential  and  Field  Arrows  results 


.  Relevant  Result* 
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Hn:  -2.493e-19 


a)  Subdomain  Settings 


b)  Boundary  Settings  c)  Electric  Potential  and  Field  Arrows 
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Electrical  Input  Impedance 


Integrate  Displacement  Current  Density  for  Input  Current 


III  =-(up(Jdy_smppn))*z 


Current  into  upper  plate 


ID  =-(down(Jdy_smppn))*z  Current  into  down  plate 


Admittance  is:  Y=1/Z 


Z  =  Impedance 


|Y|=l/V=abs(ID)+abs(IU) 


Magnitude 


For  V  =  1 


G=  real(ID)-real(IU) 
B=  imag(ID)-imag(IU) 


Real  Part 
Imaginary  Part 


Capacitance  is:  b/w 

C=  (imag(ID)-imag(IU))/(2*pi*freq) 


Y  (S) 


NEWPORT 


Piezoelectric  Stack 


Revolutionary  Research  .  .  .  Rclcront  Result* 

In-air  Admittance  Magnitude  Response 
Measured  ( _ )  and  COMSOL  Modeled  (oooo) 

Y 


Frequency  (kHz) 
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COMSOL  Fluid,  Shell  and  Stack 
Subdomain  Settings  and  Boundary  Settings  for 
Pressure  Acoustics  (acpr)  application 


a)  Subdomain  Settings  b)  Bound? 

Blue  -  Water  fluid  two  wavelength  at  3  kHz 

Yellow  -  PML  one  wavelength  Mesh  16312  Triangu|ar  Element 
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Sound  Pressure  Level 


Revolutionary  .  .  .  Releront  Results 


COMSOL  SPL  is  given  by, 


v^iu  v a  a  /  a  Aey 


Pref  =  1  M^3  f°r  Water 


2  in  denominator  indicates  pressure  is  a  peak  value, 

Therefore  Voltage  Potential  on  transducer  should  be  a  peak  value  V  =  1.414 
The  more  common  form  of  SPL  in  terms  of  rms  pressure  is  given  by, 


III  UJ  I  I  t^/ 


Transmit  voltage  response  for  2-D  model  is  given  by, 

TVR  =  SPL- 1 0log  1 0(  1/h)- 1 0log  1 0(  1/r)  dB  ref  1\xPa/Vrms  at  1m 


h  is  the  transducer  height  in  the  z-direction 
ris  the  radial  distance  pressure  point 
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TVR  (dB  re.  1uPaA/@1m) 
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Transmit  Voltage  Response 
Measured  ( _ )  and  COMSOL  Modeled  (oooo) 

28-Aug-96  11:53:04 


Transmit  Voltage  Response  Stave  Side 
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Measured  ( 

28-Aug-96  11:53:04 


In-water  Conductance  Response 

_ )  and  COMSOL  Modeled  (oooo) 


CONDUCTANCE 


B  (uS) 


WARFARE  CENTERS 


NEWPORT 


28-Aug- 


ln-water  Susceptance  Response 

( _ )  and  COMSOL  Modeled  (oooo) 


Revolutionary  .  .  .  Rclcrant  Rebuff  * 


SUSCEPTANCE 
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Beam  Patterns  in  Omni  Mode 
Measured  and  COMSOL  Modeled 
at  3  kHz,  5  kHz  and  9.6  kHz 


Research  .  .  .  Rclcrant  Result* 


SENECA  LAKE 

28-Aug-96  12:03:31 

9600  Hz  File:  NUWCDOG.16 

CML  288  TAP 

58FT  100VRMS 

F-37  #280  SS-206.6 

S/N  001  (PATTERN) 

Beam  Pattern 

OdB  ref  =  193  dB  re.  uPa  @  1m 
-3dB  BeamWidth  =  35.36  Deg. 

Dl  =  3.11  dB 


Beam  Pattern 

OdB  ref  =  110.3  dB  re.  uPa  @  1m 
-3dB  BeamWidth  =  38.63  Deg. 

Dl  =  3.2  dB 


b)  COMSOL  Modeled 
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Directional  Dogbone 
Flextensional  Sonar  Transducer 


Revolutionary  .  .  .  Releront  Results 
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FEA  Shell  and  Stack  Modes 


Revolutionary  Research  .  .  .  Relerant  Results 


Omni  Mode 


Dipole  Mode 


Directional  Mode 


Ea=1V 
Eb=1  V 


Ea=1V 

Eb=-1V 


Ea/Eb  =  Complex  Drive 


Displacements  and  PZT  Stack  Electric  Fields 
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Voltage  Drive  Configuration 


Complex  Far  Field  Pressure 


+ 


p 


dirct 


Complex 
Drive  Voltage 
Coeff 


E„=(Pd  +  P0)IPd=\  +  R 
Eb  =  (Pd-PJ/Pd=l-R 


vtfiere  R  =  PJ  P, 
E/E,=l  +  R/l-R 


All  Pressures  Normalized  to  Dipole  Pressure  Field  Pd 
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COMSOL  Modeled 


Research  .  .  .  Relevant  Result* 


a)  Pressure  and  b)  SPL  Surface  Plots  in  the 
Omnidirectional,  Dipole  and  Directional  Modes  at  3  kHz 


Omni  Mode  Dipole  Mode  Directional  Mode 


a)  Pressure  Surface  Plots 


b)  SPL  Surface  Plots 
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FEA  Molded  Beam  Pattern 

at  3  kHz 


Omni  Mode  Dipole  Mode 


Directional  Mode 


EA=1V 

EB=1V 


EA=1V 

EB=-1V 


EA=  -0.721 -j0.244= 

0.761  @-161 .3  deg 
EB=1  V 

Front/Back  ratio  =  67.6  dB 


Same  results  obtained  from  1kHz  -  6  kHz,  with  Front/Back  >  50  dB  over  Band 
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Summary 


Revolutionary  .  .  .  Releront  Results 


□  The  Directional  Dogbone  Flextensional  Transducer 

•  Generates  Cardioid  Directional  Beam  Patterns 

»  Front  to  Back  ratio  >  50  dB 

»  That  can  be  steered  Left  or  Right 

-  Reduces  Ambiguity 

-  Reduced  reverberation 

-  Improve  detection  rates 

•  Broadband  >  Octave 

•  Reduced  Size,  Weight  and  Cost 
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